Constraints on flavor-changing Z' models by Bs mixing, Z' production, and Bs —>■ jjl^ ^ 
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Certain string-inspired Z' models have non-universal interactions to three families of fermions and 
induced tree-level flavor-changing couplings. We use recent results on Bs-Bs mixing to constrain 
the size of the flavor-changing couplings in the b-s sector. In some highly predictive Z' models, 
such a constraint on b-s coupling can be translated into the flavor-diagonal couplings. Based on the 
Z' production limits at the Tevatron, we obtain the limit on the leptonic couplings of Z' and then 
make predictions for Bs -^ t'-^ ^"^ branching ratios. We conclude that with the present constraints 
from Bs mixing and Z' production, the muonic decay of Bs may not be observed at the Tevatron 
if the projected integrated luminosity is less than 0(5 — 10) fb. 



Searches for flavor-changing neutral currents (FCNC) 
have been pursued for many years. So far, the sizes 
of FCNC in the u-c, b-s, s-d, and b-d sectors, in gen- 
eral, agree with the Standard Model (SM) predictions, 
namely, those derived from the Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism in higher order. The FCNC 
effect in b-s sector was recently confirmed in the Bs me- 
son mixing observed by both CDF and D0 : 

CDF : AM, = 17.33 +"ail (stat.) ± 0.07 (syst.) ps"^ , 



D0 



AMs = 19.0 ± 1.215 ps" 



where we have converted the 90% C.L. bound 17 < 
AM, < 21 ps"'^ of D0 into la range assuming the error 
is Gaussian. We combine both results, again assuming 
Gaussian errors, and get 



AM^P = 17.46 tltl PS-' 



(Icr range) . (1) 



We re-evaluate the SM prediction ^1], using the best- 
fitted inputs (given later) before the announcement of 
the new B^ mixing data. 



AM™ = 19.52 ± 5.28 ps~^ 



(2) 



It is important to use the best-fitted inputs without the 
new Bs mixing data in order to determine if there is 
any discrepancy between the data and the SM predic- 
tion. Measurement of Bs mixing is often used to deter- 
mine the value of \Vts\, but this is clearly inappropriate 
when the mixing has additional contribution from new 
physics. The SM prediction in Eq. (j^jl contains large un- 
certainty from the hadronic parameters, nevertheless, the 
data agrees fairly well with the SM value. Therefore, we 
can use the Bs data to constrain new physics that may 
induce the b-s transitions. 

Another important channel to search for FCNC is the 
muonic Bg decay, Bs -^ M"*"/^") which has the largest 
chance to be detected at hadronic machines. In the SM, 
this process is loop-suppressed. However, many exten- 
sions of the SM predict a branching ratio large enough 



to be seen at hadron colliders. We consider an FCNC Z' 
model inspired by string theory in this letter. The Bs 
mixing and Bs -^ M"*^/^" ^r^ highly correlated because 
from the Bs data one can constrain the FCNC b-s-Z' 
coupling, which is an essential element in the calculation 
of the muonic decay. One additional element is the ^-/x- 
Z' coupling. In order to make a reliable prediction for the 
muonic decay branching ratio, we take into account the 
(t(Z') • B{Z' -^ e+e") limits from the Tevatron. In the 
Z' model considered here, the FCNC b-s-Z' coupHng is 
related to the flavor-diagonal couplings qqZ' in a predic- 
tive way, which are then used to obtain the upper limits 
on the leptonic UZ' couplings. Therefore, we are able 
to predict the maximally allowed branching ratio for the 
muonic decay of Bs- The predicted branching ratio is 
always less than 9 x 10"^ for Mz' = 200 - 900 GeV. It 
implies that with the present constraints from Bs mix- 
ing and Z' production, the muonic decay of Bs may not 
be observable at the Tevatron if the projected integrated 
luminosity is less than 0(5 — 10) fb. 

In some string-inspired models, the three generations 
of SM fermions are constructed differently and may re- 
sult in family non-universal couplings to an extra U{1) 
gauge boson, Z' . Without loss of generality, we consider 
the case that the Z' couples with a different strength to 
the third generation, as motivated by a particular class of 
string models [3| ■ Once we do a unitary rotation from the 
interaction basis to mass eigenbasis, tree-level FCNCs are 
induced naturally. Several works have recently been done 
regarding the FCNCs in the down-quark sector [li y, Ij] . 
In order to increase the predictive power, we assume that 
the left-handed (LH) up-type sector is already in diago- 
nal form, such that Vckm = VdL, where VdL is the LH 
down-type sector unitary rotation matrix. Since we do 
not have much information about either the right-handed 
(RH) up-type or the RH down-type sectors, we simply as- 
sume that their interactions with Z' are family-universal 
and flavor-diagonal in the interaction basis. In this case, 
unitary rotations keep the RH couplings flavor-diagonal. 



Therefore, the FCNCs only arise in the LH d-s-b sector. 
The coupUngs depend on the CKM matrix elements and 
one additional parameter x, which denotes the strength 
of Z' coupling to the third generation LH quarks relative 
to the first two generations. Consequently, if a; is an 0(1) 
parameter but not exactly equal to 1, the b-s-Z' coupling 
will induce a significant FCNC effect. 

We follow closely the formalism given in Ref. pj- We 
assume for simplicity that there is no mixing between Z 
and Z'^ as favored by the precision data. The current 
associated with the additional U{1) gauge symmetry is 






.(2) 



g(2) 



V'j 



(3) 



So far, we have not specified the RH chiral couplings of 
the down sector nor the chiral couplings of the up sector. 
In order to obtain constraints from Z' production at the 
Tevatron, we take the following assumptions 



lOfll^lQil, \QlR\ = \QiM\ 



(7) 



Such assumptions are reasonable; many Z' models pre- 
dict chiral couplings to be of a similar order {e.g., the Z^ 
model has all chiral couplings equal to 1/V24). More- 
over, the prediction of Z' production depends on the fac- 
tors Ql^" + Q'j^ . Thus, changing to another Z' model 
would not affect the limits significantly. 

Within the SM the mass difference in the Bs system is 



where e. 



(2) 

■>l>L,Rij 



is the chiral coupling of Z' with fermions i 



and j running over all quarks and leptons. The Z' cou- 
plings to the leptons and up-type quarks are assumed 
flavor-diagonal and family-universal: e^ ^ = Q^ ^1, 
^L B.~ Ql fll ^^'^ ^L — Ql^ where 1 is the 3x3 identity 
matrix in the generation space and Q^ ^, Qi r and Q1 
are the chiral charges. On the other hand, the interaction 
of Z' with the down-type quarks is 



4'^ = -.92^; {d, s, b)^^ {eiPL + eiPR) I s 

b 

' '(4) 

where the subscript / denotes the interaction basis. For 
definiteness in our predictions, we assume 



1 

ei = gi I 1 
X 



^R — Qb 




(5) 



The deviation from family universality and thus the mag- 
nitude of FCNC are characterized by the parameter x in 
the b^-bh-Z' entry. The chiral charges have to be speci- 
fied by the Z' model of interest. 

When diagonalizing the down-type Yukawa matrix, we 
rotate the LH and RH fields by VdL and VdR, respectively. 
With the form of e^ ^ assumed as in Eq. ©, the RH 
sector remains flavor-diagonal in the mass eigenbasis but 
the LH sector V'j^e^Vdi is in general non-diagonal. With 
no mixing in up-quark sector, we have Vckm = VdL, 
making the model very predictive. Explicitly, 



B 



L = ^dL'^L^dL - VcKM^L^CKM 



(6) 



Qi 



1 (a;- 

i)VtdVt* 
mtdV.l (x- 



1 {x- 



mtbVti 



67r2' 



AM|f = ^M^mBjUVtbV,:)'v2BSo{xt) 



x[a.(mfc)]-6/23 



An 



■h 



BsArrib) ,(8) 



where So{xt) = 2.463 and the NLO short-distance QCD 
corrections are encoded in the parameters r]2B — 0.551 
and Js ~ 1.627 |6|. We have taken Mr = 5.3696±0.0024 
GeV, and tb, = 1.466 ± 0.059 ps'^ Q. For the Wolfen- 
stein parameters, we use the CKMfitter results after EPS 
2005: A = 0.22622 ± 0.00100, and A = 0.825lJJ;^ig , 
p = 0.207to°Qlt, and fj = 0.340 ± 0.023 0. The hadronic 



parameter fB^y Bb^ = 0.262 ± 0.035 is taken from the 
lattice calculation [2j. After taking the mean for asym- 
metric errors, we obtain the SM prediction ^ as in Eq. 1^)1. 
The effect of LH FCNC induced by Z' is given by 



AMSm 



1 + 3.57X lO'^(pf)%2'0i 



= 0.894±0.243 



(9) 
where (j)f^ is the weak phase associated with the coupling 
Bf". In the model that we consider 
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92Mz ^,b 



giMz' 



-Bf 



92 Mz 



giMz' 



{x - mivtbv:, 



(10) 



and (f)"^ = 180°. We have combined the relative errors of 
AMfP and AAf™ in quadrature in Eq. ®. 

We show the allowed parameter space of {pf^,<j>f^) in 
Fig. n For (pj^ = or 180°, p|f is constrained to be less 
than 6.20 x 10^^. In more general models, (j>f^ may have 
a different value. For example, if (f>j^ — 90°, pj^ is con- 
strained to be less than 9.87 x 10~^. Note that there are 
regions with pf^ > 9.87 x 10~* also allowed by the current 



'tb 



where we have made simplifications using the unitarity 
of Vckm- It is interesting to note that the sizes of the 
flavor-changing couplings satisfy the hierarchy: ji?^"! > 
> 



K' 



\Bf\ 



^ For consistency, this SM value is obtained without referring to 
the measured AM^j because Z' can also have contributions in the 
B^ system, even though the uncertainty in Eq. |2j could have 
been much smaller by doing so. 
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FIG. 1: The allowed parameter space in the Z' model with 
FCNC only in the LH sector. The shaded area corresponds 
to the Icr C. L. limits of l\MT^ j l\Mf^ = 0.894 ±0.243. The 
central value 0.894 corresponding to t\Mf^^ = 17.46 ps~^ is 
also indicated. 



FIG. 2; The upper limits on gi\J {'Q'\} +Q%^)/2 obtained 
using the CDF 95% C.L. limits of cr(Z') ■ B{Z' ^ e+e") 
tabulated in Table |l] The constrained value of pf^ — 6.20 x 
10^* is used. 



AMs constraint, e.g., 2.15 x lO"^ < pf < 2.45 x 10"^ foj. 
(pfl = 90°. However, some of these regions correspond to 
Z' contributions larger than the SM contributions. Al- 
though not completely impossible, we think it is unlikely 
and thus leave it out from the discussions in the rest of 
the paper. 

The production cross section of Z' followed by the lep- 
tonic decay is given by 



gt iMz 



a{pP^Z'^(.+n = i^- 



144 s r^ 
dx 



Ql+QI 



-u.d.s.c 



(of + oiO/ =^'(-)aG) 



(11) 



where y/s = 1960 GeV, r ~ M^,/s and Tz' is the total 
width. The partial width Z' -^ f f is 



r{z' ^ //) 



487r 



yr^v 



2{\Qi\' + \Qi\' 



TABLE L The 95% C.L. limits on a{Z') ■ B{Z' -^ e+e") 
given by the preliminary CDF result in Ref. jlO|] as a function 
of Mz' . 



Mz' (GeV) 


a ■ 5^5 (pb) 


Mz' (GeV) 


a ■ £95 (pb) 


200 


0.0505 


600 


0.0132 


250 


0.0743 


650 


0.0136 


300 


0.0289 


700 


0.0134 


350 


0.0404 


750 


0.0126 


400 


0.0261 


800 


0.0171 


450 


0.0259 


850 


0.0172 


500 


0.0172 


900 


0.0215 


550 


0.0138 


950 


0.0246 



x(l 



fi) + 12fiQ{Q{ 



(12) 



where Nf = 3(1) for quark (lepton) and fi 



\/Nf 



We have included all leptonic and hadronic modes in the 
total width. Note that the FCNC contributions are neg- 
ligible and Z' -^ W^W" is highly suppressed by the 
Z-Z' mixing angle, which is severely constrained by elec- 
troweak precision data [ll| . 

The most recent (though preliminary) upper limits 
on the Z' search was performed by CDF with an inte- 
grated luminosity 819 pb^-'^ ,lj|. We read off the limits 
of a ■ B{Z' —> e^e~) from their figure, and tabulated in 
TableQ] We use the constrained value of p|J' to obtain the 
value for g2Q'L for each Mz', which is in turned related 
to 52 Qi by Eq. ©. With our assumptions in Eq. {TJ, we 

can then obtain the upper limits on g2\/ {Ql^ + Q'r^)/'^ 
from the CDF 95% C.L. li mits of a{Z') ■ B {Z' -^ e+e^). 
We show the Hmits of 52\/(Qi^ + Q%'^)/2 with the con- 
strained value of pf' — 6.20 x 10^* in Fig. 12 At smaller 
Mz' the limits are insensitive to the value of x, while the 
difference is more visible at larger Mz'. When x gets 
larger toward 1, the allowed value for Q^ from p^^ in- 
creases. Thus, the corresponding upper limits on Q|^ ^ 
have to be smaller. That is why at each Mz' the cho- 
sen values of a; = 0.1 — 0.9 go from top to bottom. The 
general trend for increasing Mz' is the increase in the 
upper limit of Q\ ^. This is easy to understand because 
at large Mz' the dominant factor in the production cross 
section is the parton density, which becomes very small 
at large momentum fractions. 

We are ready to compute the maximally allowed decay 
rate for Bg -^ IJ'~^I^~ ■ We ignore the RG running effect at 
the b-s-Z' vertex, which is good enough for an order-of- 
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FIG. 3: The allowed parameter space on the {pt,p1^ = p^**) 
plane using the current upper bounds on B{Bs — > pL^ p~): 
1.0 X IQ-'' (CDF) and 2.3 x 10"'' (D0). The upper limit of 
pf^ for (j)"^ — 0° and 90° are also indicated using the blue 
dashed and red dashed lines, respectively. 



magnitude estimate. The branching ratio of Bs — > /i^A* 
is given by 
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FIG. 4: Dependence of the upper limits of B{Bs —^ p^ p, ) 
on Mz' , where different sets of lines correspond to different 



values of x varying from 0.1 to 0.9. We have set p^'' 



PZ 



\VtlVts\ 



(13) 



where 



UM _ 92 Mz ^e 
Pl,r - „ ,,r '^L,R 



(14) 



at the weak scale. One can find the definition of 
Y{m^/M'^) in the SM part in Ref. 0; its value is 
about 1.05 here. Using the central value of the av- 
eraged Bs lifetime tb^ = 1.461 ps, corresponding to 
Fb, = (4.49 ± 0.18) X 10-13 GeV and /b, = 230 MeV, 
we obtain a SM branching ratio of about 4.2 x 10"^. The 
current upper limits on B{Bg —> jX^ pT) from CDF and 
D0 based on 780 and 700 pb"-^ data, respectively, are 
El: 



values of the parameter x in our model, ranging from 0.1 
to 0.9. It is seen again that there is less variations in 
the predictions for small Mz' ■ The result that the maxi- 
mally allowed branching ratio is always less than 9 x 10~^ 
implies that Tevatron may not have the capability to ob- 
serve the muonic decay of B^. These upper limits are 
at least one order of magnitude below the current best 
upper bound on B(Bs —> pL~^ p,^) quoted above. 

We conclude that with the present constraints from 
Bs mixing and Z' production, the muonic decay of Bs 
may not be observed at the Tevatron if the projected 
integrated luminosity is less than 0(5 — 10) fb. However, 
at LHCb, with anticipated production of 10^^ bb per year, 
the expected branching ratio of order 10~® is observable. 

This research was supported in part by the National 
Science Council of Taiwan R.O.C. under Grant Nos. NSC 
94-2112-M-007-010- and NSC 94-2112-M-008-023-, and 
by the National Center for Theoretical Sciences, and in 
part by the U.S. Department of Energy under Grants No 
DE-FG02-96ER40969. 



B{Bs^p+p^) < 1.0 xlQ-^ (CDF) 
B{Bs -^ p+p-) < 2.3 X 10-^ (D0 ) . 

Given the assumption that RH FCNC couplings are ig- 
nored and the fact that (j)^^ — 180° in the current model, 
we plot the allowed parameter space on the pf^-p^^ plane 
in Fig. 13 We have taken p'^ = p'^'^. 

If we take the upper limit of pf^ — 6.20 x 10"'', and 
the conservative upper limit of g2Qi given in Fig. [3 for 
different values of Mz' , we obtain the upper limits of the 
Bg — > p^ p~ branching ratio, as shown in Fig. 0] Dif- 
ferent sets of lines in the drawing correspond to different 
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